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ABSTRACT 

We present Doppler images of the young K5V - K7V rapid rotator LO Peg from seven 
nights of continuous spectroscopy obtained in 1998 from July 04 to July 10. The images 
reveal the presence of a strong polar cap with appendages extending to mid-latitudes, but no 
starspots are seen below 15°. We briefly discuss the distribution of spots in light of recent flux 
transport simulations which are able to reproduce the observed latitude dependence. 

With the full timeseries of spectra, of which 314 are useful, many phases are observed 
three times over the seven nights of observations. Using starspots as tracers of a solar-like 
latitudinal differential rotation in our image reconstructions, we find that the equatorial regions 
complete one more rotation than the polar regions every 181 ±35 d. LO Peg is the second 
coolest star for which such a measurement has been made using indirect imaging methods. 
The degree of latitudinal shear is less than that seen in G and early K dwarfs, suggesting a 
trend in which differential rotation decreases with stellar mass in (pre-)main sequence objects. 

Key words: Line: profiles - Methods: data analysis - Techniques: miscellaneous - Star: LO 
Peg (BD +22 4409) - Stars: activity - Stars: atmospheres - Stars: late-type 



1 INTRODUCTION 

Several methods which enable us to measure the surface differen- 
tial rotation in late-type stars with convective envelopes have ap- 
peared in the literature. Variability of stellar activity diagnostics 
such as the Ca II H & K emission cores has been used to deter- 
mine stellar rotation periods as exemplified by the long term Mount 
Wilson Survey |WilsQril 978; Baliuna s"et alll995l) . In many stars, 
the measured rotation period is found to vary periodically from 
one epoch to another over timescales similar to the 1 1 year solar 
magnetic activity cycle. If we adopt the solar paradigm, in which 
changes in latitude of the active regions occur throughout the cycle, 
we can attribute the measured period changes to a latitude depen- 
dent rotation in stars. The scatter in measured periods gives some 
indication of the degree of differential rotation shear on the star 
iDonahue. Saar & BaliunaJl99fj) . 

A similar method based on photometric modulation due to 
starspots, in a sample of stars, has been foun d to yield a trend in the 
degre e of surface shear with rotation rate iHalll99ltlHenrv et alJ 
1995). Long term variations of the photometric period been shown 
to be indicative of m agnetic activity cycles (Messina & Guinarl 
120021) . More recently, iReiners & Schmittl <2003l) have found that 
main-sequence stars of spectral types F and G exhibit high degrees 
of differential rotation. The rotational shear manifests itself in the 
Fourier transforms of the rotationally broadened line profiles. This 
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allows precise measurements of shear from line profiles alone for 
stars with intermediate rotation rates. 

A further method of determining the latitude dependence of 
rotation is to use spot distributions from indirect imaging tech- 
niques. Here we do not have to make assumptions about the lat- 
itude extent of starspots in order to determine the magnitude of 
shear between the equator and pole. We are however restricted 
to observing rapid rotators if we are to obtain sufficient reso- 
lution elements across a rotationally broadened profile to make 
this technique viable. The number of resolution elements is deter- 
mined by a number of stellar broadening mechanisms, but is ulti- 
mately limited by the instrumental resolution. By cross-correlating 
constant latitude strips in Doppler images derived at closely re- 
lated e pochs, the latitude dependent sur face shear can be deter- 
mined iDonati & Collier Cameronlll997t) . Alternatively, a differ- 
ential rotation law (usually assumed to be solar-like) can be in- 
corporated into the imaging model, and the whole data set from 
several nights combined to derive an image iDonati et alj|2000l 
IPetit. Donati & Collier Cameronl2002T) . 

LO Peg (BD+22 4409) is among the least massive of the 
young rapid rotators in the solar neighbourhood. It was detected by 
the ROSAT WFC EUV all-sky survey as the source RE J2131+23, 
and by the Extreme Ultraviolet Explo r er survey as the source 
EUVE J2131+23.3 iMalina et alJll994l) . Ijeffries & Jewed <1993l) 
identified LO Peg as a member of the Local Association on the 
basis of its Galactic space motions and a large EUV to bolometric 
flux ratio. LO Peg was first studied in detail by Je ffries et alJ<1994 
who determined an axial rotation period of 0.42375 d from V-band 
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Table 1. Journal of observations of LO Peg and standards on 1998 July 04 - 10 



Object 


UT Start 


UTEnd 


Exp time [s] 


No. of frames 


Comments 








1998 Jul 04 






Gl 673 


22:27 




100 


3 


K5V template 


LOPeg 


01:32 


05:57 


200 


49 










1998 Jul 05 






Moon 


21:57 


22:07 


200, 30, 60 


3 


Solar template 


LOPeg 


01:46 


05:56 


200 


47 










1998 Jul 06 






Gl 687 


21:38 


21:53 


300, 200, 200 


3 


M3V template 


Gl 649 


22:03 


22:17 


200 


3 


Ml V template 


LOPeg 


01:40 


05:54 


200 


47 










1998 Jul 07 






LOPeg 


01:46 


06:01 


200 


48 










1998 Jul 08 






LOPeg 


01:34 


05:59 


200 


50 










1998 Jul 09 






LOPeg 


01:51 


06:00 


200 


47 










1998 Jul 10 






G1447 


20:49 


21:23 


600 


3 


M4V template 


LOPeg 


01:41 


06:00 


200 


47 





photometric observations. It was also noted that a slightly shorter 
period of 0.38417 d gave a physically acceptable solution with low 
false alarm probability and folded lightcurve with little scatter. 
Photometry yielded a visual magnitude of 1/ = 9.19±0.05 
and colours (B — V) = 1.08 ± 0.02, (U — B) = 0.82 ± 0.08, 
{V - R)kc = 0.59 ± 0.02 and (V — I)kc = 1-23 ± 0.03. 
Study of the Lithium 6708 A line yielded an abundance of N(Li) 
= 1.30 ± 0.25, while the iron abundance [Fe/H] = 0.0 ± 0.1. 
The equatorial rotati on velocity wa s determined to be v sin i = 
69 ± 1 tans^. ljeffriesetalJ 1 1994) also obtained a model fit to 
the spectrum of LO Peg, and found a spectral type of K5 - K7 with 
an age in excess of 30 Myr. 

I Jeffries et alJ J 1994) also determined the axial inclination of 
LO Peg to be 50°. By studying the Ha line, which appears in 
emission, these authors showed that LO Peg does not exhibit 
prominence features, which appear as Ha absorption transients 
in other si milar stars as the cool optically thick clouds cross the 
stellar disc ICollier Cameron & Robinsonll 19891) . The lack of ob- 
servational evidence indicates that if prominences do exist, they 
must be present only above low latitudes where they will not 
cross the s t ellar d isc in the observer's line of sight. More recently, 
lEibe et all 1 19991) found that variability of the Ha line indicates 
str ong d ownflows of absor bing material. A p revious imaging study 
by iLister. Collier Cameron & Bartusl ll999l) has yielded surface 
Doppler images showing a polar cap and a low latitude band of 
features on the surface. Observations at two epochs allowed these 
authors to distinguish between the two proposed axial rotation pe- 
riods discussed above, indicating that the longer 0.42375 d period 
is preferred. 

The K5V -K7V spectral type make LO Peg an important ob- 
ject because no other single stars of this mass have been studied. 
In fact, no stars between K3V and and M1V have been imaged 
to date. At this spectral type, the radiative core decreases in size 



as the star approaches the fully convective regime at spectral type 
early to mid-M. If a solar-like interface dynamo is at work, accord- 
ing to Schiissler et al. 1 1996) and Granzer et al. 1 2000) we may ex- 
pect only intermediate to high latitude eruption of magnetic flux. 
In this paper, we present images derived from data obtained over 
seven nights and determine the differential rotation for LO Peg. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations were made with the Multi Site Continuous 
Spectroscopy (MUSICO S) instrument (resolution ~ 30000) 
(Baudrand & Bohm 1992) at the Isaac Newton Telescope on 1998 
July 04 - 10 (TableG). The seeing conditions were generally good 
at ~ 1" with attenuation due to atmospheric dust on the nights of 
July 07 & 08. The 2Kx2K SITel CCD was used in conjunction 
with the Cassegrain fibre fed spectrograph, giving a useful wave- 
length coverage of 4359 A to 7352 A in which a total of 51 orders 
were recorded. 

Pixel to pixel variations were removed using flat-field expo- 
sures taken with an internal tungsten reference lamp. The worst 
cosmic ray events were removed at the pre -extraction stage using 
the FIGARO routine BCLEAN JShortridgelll993l) . Scattered light 
was modelled by fitting polynomials of degree 7 to the sets of 
inter-order pixels at each X-position in each frame. The spectra 
were extracted usin g ECH OMOP, the echelle reduction package 
developed by Mills 1 1994). The Thorium- Argon arc-frames used 
for wavelength calibration were extracted in conjunction with a 
target spectrum, and calibrated using this package. The orders 
were extracted using ECHOMOP's implementation of the optimal 
extraction algorithm developed by Home] Jl98rJ) . ECHOMOP 
propagates error information based on photon statistics and readout 
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Figure 1. Time series spectra for 1998 July 04- 10, with phase plotted against velocity. The mean profile has been subtracted from each time series. White 
features correspond to starspot signatures. 



noise throughout the extraction process. 2.1 Least squares deconvolution 

We apply least squares deconvolution iDonati et all 19971) in order 
to obtain sufficient S/N in the rotation profile at each phase, en- 
abling us to reliably and consistently reconstruct surface images. 
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Figure 2. Maximum entropy regularised reconstruction of LO Peg for (a) 1998 July 04 & 05, (b) July 06 & 07 and (c) July 08, 09 & 10. 

Due to extreme vignetting in the faintest (blue) orders, we applied of typically 750 for each observation. We used a velocity bin size 

a mask to exclude the regions where few counts were recorded. of 4.5 km s _1 chosen to match the CCD pixel size at the mean 

Nevertheless, 3897 images of 2513 lines (due to overlapping wave- weighted central wavelength of the deconvolved line (5399 A), 
length ranges in adjacent orders) were used to derive single high 
S/N ratio profiles. A mean input S/N ratio of 15 was boosted by a 
multiplex gain of 50, yielding mean absorption profiles with S/N 
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Figure 3. Maximum entropy regularised reconstruction of LO Peg for (a) 1993 August 04 & 05, (b) 1998 all nights combined without differential rotation and 
(c) 1998 all nights combined with differential rotation. 



2.2 System parameters in light of differences with the values found for these parameters by 

other authors. 

Stellar system parameters were determined empirically from the 
full data set by searching for the combination of parameters which 
minimised x 2 (Table|2j- We discuss the results further in Section lT2l 
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Table 2. System parameters for LO Peg. 



P[d] 

v r [kms - 1 ] 
v sin i [kms ] 
Axial inclination [deg] 



0.423229 ± 0.000048 
-20.90 ± 0.25 
65.84 ± 0.06 
45.0 ± 2.5 



3 RESULTS 
3.1 The Images 

The maps of LO Peg from 1998 July were prod uced with our im- 
age c ode, DoTS (Doppler Tomography of Stars) ICollier Cameronl 
1200 ll) and are presented in Figure|2| Since the phase coverage in a 
single night was approximately 0.4, we have combined either two 
or three nights in order to obtain images with good phase cover- 
age. This allows us to compare image consistency and to check 
for surface evolution of starspots. Immediately obvious in all the 
images is the very strong polar cap which shows reproducible mor- 
phology from one night to the next. Since the axial inclination of 
LO Peg is 45°, stellar latitudes above this value are visible at all 
times. In essence this ensures reasonable reliability of reconstructed 
starspots above 45° for phases where no observations were made, 
although projection effects and limb darkening may mean that the 
morphology is less well constrained than for well observed phases. 
Although the polar cap reaches a radius of 25° at some longitudes, 
it is by no means uniform. In this respect, there is evidence for 
lower spot filling, very close to the pole in the July 06 & 07 and 
July 08, 09 & 10 images. The high latitude spot filling may possi- 
bly be due to two large spots or spot groups located approximately 
180° apart, and centred around phases 0.0 to 0.1 and 0.5 to 0.6. 

Below latitude 60°, a number of features are repeated in all 
the images. Intriguingly, most of the starspot structure at this time 
appears at, or above 45° with only a few small spots or spot groups 
below this latitude. In fact, no spots are seen below 15° at the time 
of observations. The mean spot latitudinal filling plot on the right 
of each image shows no real evidence for preferred latitudes below 
the obvious polar feature as is often seen on other stars of earlier 
spectral type (e.g. B arnes et all200ll) . 



3.2 Comparison with 1993 images of LO Peg 

Our im ages show similarities with those presented bv lLister et alJ 
( 1999) while at the same time revealing marked differences. The 
polar cap is somewhat weaker, but shows the same non-uniformity 
as in the 1993 images. The most notable difference is in the location 
of the low latitude features in 1993 which are confined to a definite 
band centred around 20° . 

Caution is neede d however when comparing images from the 
two epochs. iLister etalJ 119991) derived a significantly different 
v sin i for LO Peg, of 69 km s~ in agreement with the measure- 
ment made by Jeffries et alJ ll994l) . Lister et al. initially estimated 
v sin i by minimising \ 2 with respect to v sin i for an unspotted 
stellar image. They also performed reconstructions (to a fixed tar- 
get x 2 ) at a range of v sin i values, and selected the value that min- 
imised the total area of spots in the image. Use of initial \ 2 can 
only be expected to give an approximate fit to the data because an 
unspotted star yields model profiles which neglect the effects of 
starspots. Starspots can alter the shape of the profile, leading to in- 
correct determination of parameters such as radial velocity, equiv- 
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Figure 4. The mean residual plots from fits to the 1998 July 04 -10 data 
using vsini = 65.84 kms - 1 and v&mi = 69 kms -1 . The dashed lines 
indicate v sin i = 65 . 84 km s — 1 . 



alent width and v sin i, if their effects are not considered. Radial 
velocity is the least sensitive parameter in this respect since the 
mean profile of a timeseries covering a complete rotation cycle will 
be symmetric. By contrast, when a strong polar spot is present, as is 
found on LO Peg, the rotationally broadened profile becomes either 
flat bottomed or more rounded, with a reduced equivalent width. 
By fitting the profiles derived from an unspotted stellar image, we 
could expect to find an optimum fit where the line equivalent width 
is too low and the v sin i too high. 

Fitting the model to the data for a fixed level of \ 2 and mea- 
suring the starspot area (a proxy for the image entropy), may also 
lead to incorrect choice of system parameters. If a fixed x 2 is cho- 
sen, the combinations of parameters yielding the minimum spot 
area may occur before many distortions in the spectra have been 
properly fitted, leading to a non-optimal choice of parameters. 

We therefore determine the system parameters by minimis- 
ing the value of \ 2 attained after a fixed number of maximum 
entropy regularised iterations. The results are presented in Table 
|2|(we do not tabulate the equivalent width which depends on the 
strength and number of lines used in the least squares deconvolu- 
tion). Si nce we do not reach a r educed \ 2 - L we use the TEST 
statistic I Skilling & Brvan 1984) to determine the value of \ 2 at 
which we begin fitting noise. Reasons for not attaining \ 2 = 1 can 
be attributed to incorrectly determined error bars (i.e. systemati- 
cally too large or too small), and also to non-flat continuum on 
either side of the deconvolved rotation profiles. While the latter 
has negligible effect on the reconstructed data, incorrectly deter- 
mined error bars have no effect other than to alter the optimal \ 2 
which can be achieved. 1-ct uncertainties for each parameter are 
determined for A^ 2 = 1 + d/Xmim where d is the number of fit- 
ted data points and x 2 min the lowest achieved \ 2 ■ This assumes 
that mis-determination of noise is the only reason that ^ 2 = 1 is not 
achieved and will thus yield conservative estimates for parameter 
uncerta inties. 

In lBarnes et alJ 1200 ll) w e also re-determine d system param- 
eters from the original study iBarnes et alJll998T) of the G dwarf 
He 699 using the same method as in the present study. There, 
we found that using spot area as a proxy for the image en- 
tropy tended to overestimate the v sin i and equivalent width, 
a finding shown in Barnes 1 1999). In the case of He 699, 
a v sini = 96 km s -1 was found using the spot area criterion, 
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compared with v sin i = 93.5 km s" 1 using the \ 2 criterion. At 
65.84 km s , the present result for LO Peg is 4.6% lower, hav- 
ing important effects on the reconstructed image. Overestimation of 
v sin i means that the radial acceleration of starspot features cannot 
be correctly determined and has resulted in starspots being placed 
in a low latitude band. 

Table |2| shows the optimal axial inclination to be 45 ° . This 
is slightly lo wer than the 50° e stimated by Lleffries et alJ ll994h 
(and used by llister etal] 1 19991) ) but well within their 10° un- 
certainty. Attempts to reconstruct images at higher inclinations 
fail to fit the flat-bottomed profile correctly, because the polar re- 
gion is too foreshortened to yield the necessary flux deficit, even 
when the spot occupancy saturates at 100%. Fig.|4|is a plot of the 
mean residual line profile for maximum entropy fits using v sin i 
= 69 km s~ and 65.84 km s _1 indicating the lower residuals in 
the wings using the latter value. These wing residuals are also 
smaller for vsini = 65.84 km s _1 when a similar plot is made 
for the 1993 data set. The goodness of fit obtained for the 1998 
data set with i;sini = 69 km s _ and vsini = 65.84 km s _1 is 
X 2 min =2.82 and x^min =2.34 respectively. If we assume that 
usini = 65.84 km s _1 is the optimally determined projected rota- 
tion velocity, we note that usini = 69 km s~ falls outwith the 
99.99% confidence interval. 

We have therefore recalculated an image of the 1993 data set 
using both nights of data (Fig. [3} along with the combined image 
for the 1998 data in order to allow direct comparison of consistently 
determined images. The images from both epochs are now very 
similar in terms of the polar cap strength and location of starspots 
at lower latitudes, as indicated by the mean spot filling plot on the 
right of each map. There is marginally less mean filling in the polar 
regions from 1993, but the low latitude features in both cases only 
occur at or above ~25°. The degree of similarity of the images 
at both epochs is most likely coincidental as the chosen ephemeris 
is different for each image reconstruction. Also, the error in the 
determined period results in an uncertainty of the time of zero phase 
by up to ±5 hrs in the five year interval between observations. 



Differential rotation x landscape 



3.3 Latitude dependent rotation 

The seven night timebase of the observations enable us to investi- 
gate the evolution of surface features on LO Peg. Any given phase 
is observed approximately three times on alternate nights which 
should be sufficient to track starspot features over the timespan. 
Additionally, the images in Fig. [^indicate that most spots are suf- 
ficiently stable during the seven nights of observations to warrant 
such an investigation. We have discussed the inclusion of differen- 
tial rotation in our Doppler imaging co de in previous publications, 
a method first applied successfully by iDonati et alJ fcOOOh to the 
early G pre-main sequence star RX J1508. 6-4423. A solar type dif- 
ferential rotation law involving a shift term applied to the angular 
velocity Qq, and a latitudinal dependent shear term, AQ, are fitted 
according to Q(8) — Qo — Af2sin 2 (#). We minimise the fit in the 
same way as for other system parameters, by obtaining the lowest 
X 2 after a fixed number of iterations. 

For LO Peg, we find that fi = 14.86 ±0.0027 rad.day -1 and 
AQ = 0.034714 ±0.006692 rad.day -1 . With an axial rotation pe- 
riod, P = 0.423229 d, this implies an equivalent phase shear ampli- 
tude A<j> = 0.005525 ± 0.001065 day -1 , or an equator-pole lap 
time of 181 ± 35 d. 




Figure 5. x landscape for July 04 - 10 fitted data set. Shown are the x 
minimum, marked by a cross, and from centre outwards, the 1 -parameter 
68.3%, 2-parameter 68.3% and 99.99% confidence intervals. The dashed 
lines mark the 1-parameter 1-cr or 68.3% values for each paremeter. 



4 DISCUSSION 

4.1 Starspot distributions 

With a spectral type of K5 - K7, LO Peg is one of the latest spec- 
tral type dwarfs which have been Doppler imaged. Additionally, 
it is the only main sequence star in the mid-late K spectral type 
range to have been studied in this way. The evolution of starspots on 
the one-week timescale of the current observations is clearly quite 
small. Over the longer five year time base, the starspot morphology 
has also changed very little in terms of latitudinal distribution. The 
polar cap appears to have remained both stable and strong at both 
epochs of observation. This is intriguing because the slightly earlier 
spectral type (K3V), but more rapidly rot ating (P = 0.38 d) Speedy 
Mic showed no evidence for a polar cap iBarnes et all200ll) when 
observed at a single epoch in 1998. The early M dwarf stars, EK 
Dra (RE 1816+541) and HK Aqr have similar periods to LO Peg but 
do not possess polar caps, at more than one epoch in the latter case 
iBarnes. James & Cameronl2004l). There is however evidence that 
the young KOV star AB Dor jKurster. Schmitt & Cutispotoll 19941) 
did not possess a polar cap in 1989, in which case these large spots 
may disappear periodically. If such behaviour is (quasi-)cyclic, it 
may be related to dynamo activity, although there is as yet no con- 
vincing evidence for this. 

By simulating a s tar ten times more active than the Sun, 
ISchriiver & Title! 1200 ll) found that the higher rate of flux emer- 
gence, at around 30 times solar, resulted in the development of 
polar caps. But magnetic flux, as traced by starspots on LO Peg, 
clearly does not emerge at low latitudes in the way we see on the 
Sun. With the flux emergence patterns seen on young rapid rotators 
in mind, Mackav et al. (2004) have investigated the role of merid- 
ional flows for stars where flux emerges up to latitude 70°. In order 
to reproduce the intermingled polarity pattern as seen on the young 
rapid rotator, AB Dor, strong meridional flows of around 100 m s~ 
are required. 
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Unfortunately we do not have magnetic maps of LO Peg 
which would allow us to compare our results more closely, but the 
the flux distribution patterns do demand comparison with the simu- 
lation maps. Ma ckav et alJ l2004) have shown that much of the flux 
piles up at high latitudes, with less at intermediate latitudes. In the 
case where the rate of flux emergence shows cyclic behaviour (with 
an 1 1 year cycle), but no change in the latitude at which flux erupts 
over that cycle, fewer spots are seen at low latitudes than interme- 
diate and high latitudes. While the amount of flux at low latitudes 
varies throughout the cycle, features within a few degrees of the 
equator are present at all times. However when a solar-like variation 
of flux emergence with latitude, varying between ±70° (rather than 
the solar ±40°) to ± 5°, is simulated, there are times when no flux 
is present in the equatorial regions. A distribution resembling this 
situation is certainly seen on the present images of LO Peg where 
we see no spots below 15°. If observation can be reconciled with 
this picture of activity, long term Doppler imaging pr ojects are re- 
quired. It is worth noting however that while Baliunas et al. 1 1995) 
found period variability, indicating the presence of active regions 
at different latitudes and differential rotation, the younger stars in 
the Mount Wilson Survey tended to show no preferred period for 
magnetic cycles. 

A solar type interface dynamo predicts only spots a t in- 
termediate latitude in rapid rotators ISchiissler et alj fl996l and 
iGranzer et afl EoOO). While spots do appear at intermediate lati- 
tudes in both 1993 and 1998, there is some structure at lower lat- 
itudes, while the images are dominated by polar or circumpolar 
spots. In many cases however, in addition to polar and circumpolar 
structure, spots my be seen in bands at low latitudes. Th is is the case 
for th e young post-T Tauri star RX J1508.6-4423 iDonati et all 
2000) and Alpha Per G dwarfs iBarnes et al. 2001) for example. 
The role of meridional flows have been touched upon above and 
are known, in the solar case at least, to possess m aximum ampli- 
tude at equatorial latitudes (e.g. Ko mm. Howard & Harvevlll993h . 
Flow rates following t he solar pattern , but with greater amplitude, 
as modelled bv lMackav et alJ 12004). will then result in a lower 
chance of seeing starspots at low latitudes. Another likely reason 
for no spots being observed at low latitudes is related to the low ax- 
ial inclination of LO Peg. At 45°, spots in the equatorial latitudes 
of the star will appear severely foreshortened. If the S/N of the data 
is not sufficiently high, and only small spots exist, our imaging pro- 
gram may not find sufficient evidence for their presence. 



4.2 Differential rotation 

The differential rotation rate which we have measured for LO Peg 
is interesting in the context of other measurements. The value 
of equator-pole lap time, S = 181 ± 35 d, lies on a trend (al- 
beit sparsely populated) with spectral type. iDonati et alJ 12000) 
found S = 50 ± 10 d for t he young G2 pre-m ain-sequence star 
RXJ 1508.6-4423, while [Marsden et al. (2004) found S = 66 ± 
14 d for the IC2602 G2 star, R 58. Over a period of 8 years, 
ICollier Cameron & Donati (2002) found that the K0 dwarf AB Dor 
exhibits systematic changes in differential rotation rate from one 
season to the next, with S = 70 d to 140 d. At approximately the 
same spectral type, the near-one-day period PZ Tel iBarnes et alJ 
2000) exhibited an equator-pole lap time of 80 ± 15 d. The latest 
spectral ty pe star for which we have a differential rotation mea- 
surement iBarnes et all2004T) is the dMl-2 object, HK Aqr which 
exhibited an equator-pole lap time formally consistent with zero 
(i.e. between -1449 d and +448 d.) Although the sense of the dif- 



ferential rotation was not clear, the magnitude is essentially very 
small. 

LO Peg clearly fits into this trend of generally decreasing dif- 
ferential rotation as we descend the main sequence from early G 
through to early M. The changing AB Dor differential rotation how- 
ever shows that there can be significant change, possibly dependent 
on the phase of a magnetic cycle. In all cases, the magnitude of dif- 
ferential rotation appears, to first order at least, to be independent 
of rotation rate. The Sun for instance has an axial rotation period 
of 27 d and exhibits an equator-pole lap time of 120 d, whereas this 
value falls only by a factor of two (twice the rate of differential 
rotation) for rapidly rotating objects with l/50th the solar period. 
These two predictions of weak rotation dependence and a decrease 
with rotation rate are predicte d by the fluid circulation models of 
Kitchatinov & Riidiser 1 1999), although the exact magnitude of the 
predicted differential rotation does not match the observation. Ad- 
ditionally, for a 1 d axial rotation period K5V star, these models 
predict a poleward meridional flow rate o f ~ 10 m s" 1 , an or der 
of magnitude lower than those required by Mac kav et all 12004) in 
their simulations. Clearly it would be desirable to determine an es- 
timate of the magnitude and direction of flows observationally. 
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